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ABSTRACT 

Asthma has a strong genetic component. The final disease phenotype results from complex interactions between environment 
and multiple genes of small-to-modest effects. We investigated whether the polymorphism in genes encoding inflammatory 
mediators and cytokines is important for solving the onset and progression of asthma. We investigated whether 31 single 
nucleotide polymorphisms (SNPs) in genes encoding cytokines or monokines (interleukin [IL]-5R, matrix metalloproteinase 
[A/IMP] 8, beta! adrenergic receptor, cytotoxic T -lymphocyte-associated antigen 4, IL-3, C-reactive protein, cytochrome P450 
(CYP) 2C9, CYP3A4, a disintegrin and metalloproteinase [ADAM] 33, cysteinyl leukotriene receptor [CysLTR] 1, CysLTRl, 
eosinophilic cationic protein, glucocorticoid receptor, and leukotriene A 4 hydrolase) are related to asthma development in 206 
Japanese bronchial asthma patients and 127 healthy controls. Using multifactor dimensionality reduction (MDR), we identified 
rsl7099451 in MMP8, using a single locus model, with a mean cross-validation of 87.0%. Using a two-locus model, 
combinations of MMP8 and rs44707 in ADAM33, and MMP8 and rs40401 in IL-3, were identified, with mean cross- 
validation consistencies reaching 45.0%. Of the SNPs selected by the MDR method, rsl7099451 in MMP8 and rs40401 in IL-3 
were regarded as the most significant results in a 2 X 2 dominant model analysis. The finding that an MMP8 allele was most 
strongly related to asthma development indicates that metalloproteinase function is crucial to the airflow limitation process 
involved in this disease. 

(Allergy Rhinol 4:el32-el39, 2013; doi: 10.2500/ar.2013.4.0063) 



The main pathogenic mechanism of bronchial 
asthma involves a particular form of chronic in- 
flammation, in which the increase in eosinophils, lym- 
phocytes, and mast cells in the respiratory tract mem- 
brane and exfoliation of the respiratory epithelium 
cause bronchial hyperresponsiveness, followed by air- 
way remodeling. 1 These pathological changes result 
from the interactions between the inflammatory medi- 
ators and the various cytokines produced by the infil- 
trated cells or the structural cells. 2 

Asthma has a strong genetic component and com- 
plex interactions between the environment and multi- 
ple genes. It should therefore be possible to specify the 
risk and possibly prevent the development of the dis- 
ease by investigating genetic polymorphisms; thus, the 
investigation of polymorphisms in genes encoding in- 
flammatory mediators and cytokines is important. For 
example, genes encoding the high-affinity IgE receptor 
(Fc epsilon RI), an atopy-related gene 3 ; eosinophil cat- 
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ionic protein, an inflammation-related gene 4 ; interleu- 
kin (IL)-4 and IL-4 receptor a-chain 5 ; IL-5 receptor 
a-subunit 6 and IL-10 7 as Th2-related cytokines; and cys- 
teinyl leukotriene receptor (CysLTR) 1-related gene 8 
have been reported to be involved in asthma. Moreover, 
an a disintegrin and metalloproteinase (ADAM) 33 mu- 
tant has been reported to be related to airway responsive- 
ness and bronchial remodeling. 9 

In the study, we investigated the role played by 31 
single nucleotide polymorphisms (SNPs), located in 
genes encoding 14 cytokines or monokines, in asthma 
development and the onset of bronchial asthma in 
Japanese bronchial asthma patients. 

MATERIALS AND METHODS 
Subjects 

The study protocol was approved by the Ethics 
Committee of Fukuoka National Hospital; all partic- 
ipants received verbal and written study information 
before giving their informed consent (approval num- 
ber 23-14). 

We studied 206 Japanese subjects (102 men and 104 
women) with newly diagnosed bronchial asthma at the 
Fukuoka National Hospital and 127 controls (46 men 
and 81 women) from the Japanese general population 
(Table 1). 



e132 



Fall 2013, Vol. 4, No. 3 



Table 1 Clinical profile of study subjects 


Clinical Profiles 


Asthmatic Patients 


Normal Controls 


Number of subjects 


206 


127* 


Age (mean ± SEM; range) 


42.8 ± 1.1 (16-75) 


35.5 ± 1.0 (20-57) 


Sex (male /female) 


102/104 


46/81 


*One sample is not included in the multifactor dimensionality reduction and association analysis because of missing genotype 


ofrs!130864. 







All asthma patients were diagnosed according to the 
Global Initiative for Asthma standard with asthma- 
related symptoms (recurrent cough, wheezing, or dys- 
pnea), with bronchial hyperresponsiveness (PC20 to ace- 
tylcholine was <8 mg/ mL), and with airway reversibility 
(forced expiratory volume at 1 second (FEV^ increased 
with 12% and 200 mL or more after salbutamol inhala- 
tion). None of these individuals were current smokers, 
had chest infections, had used antibiotic or oral predni- 
sone in the previous 4 weeks, or had heart or lung disease 
other than asthma. Controls had no current disease or 
abnormal chest radiograms and were self-reported as 
healthy. 

Selection of Polymorphisms 

Using public databases, e.g., Pub Med and Online 
Mendelian Inheritance in Man, 10 we selected 14 genes 
that have been implicated in the onset of asthma or the 
asthma phenotype (Table 2), viz., IL-5R, matrix metallo- 
proteinase (MMP) 8, f32 adrenergic receptor, cytotoxic 
T-lymphocyte-associafed antigen 4, IL-3, C-reactive pro- 
tein, cytochrome P450 (CYP)2C9, CYP3A4, ADAM33, 
CysLTRl, CysLTR2, eosinophil cationic protein, gluco- 
corticoid receptor, and leukotriene A 4 hydrolase. 

Genotyping 

Venous blood (7 mL) was collected from each subject 
into tubes containing 50 mmol/L of EDTA; from these 
samples, genomic DNA was isolated with a FlexiGene 
DNA Kit (Qiagen, Tokyo, Japan). Genotypes at 31 
polymorphisms in the 14 genes were determined 
with a fluorescence or colorimetry-based allele-specific 
DNA-primer-probe assay system (Toyobo Gene Anal- 
ysis, Tsuruga, Japan). 

Polymorphic regions of each gene were determined 
using the intercalator-mediated fluorescence resonance 
energy transfer probe method (Toyobo Gene Analysis). 
In brief, the reaction was performed in a volume of 25 
fxL containing 20 ng of genomic DNA, enzyme reaction 
buffer, 2.5 mmol/L of MgCl 2 , 0.2 mmol/L of deoxy- 
nucleotide triphosphates, and 1.25 U of rTaq, contain- 
ing anti-Taqplus High (Toyobo Gene Analysis). The 
amplification protocol was performed in a PE 9700 
thermocycler (Applied Biosystems, Tokyo, Japan) and 
included an initial denaturation at 95°C for 5 minutes; 



40 or 45 cycles of denaturation at 95°C for 30 seconds, 
annealing at 65°C for 30 seconds, and primer extension 
at 72°C for 30 seconds; and a final extension at 72°C for 
2 minutes. The polymerase chain reaction products 
were then mixed with 5 mL of a solution containing 40 
mmol/L of EDTA, 20 pmol of Texas Red-labeled 
probe, and 1/2000 SYBR Green I (Invitrogen, Carlsbad, 
CA); these samples were then transferred to an ABI 
Prism 7900 instrument (Applied Biosystems, Foster 
City, CA) for measurement of the melting point as 
temperatures were increased to 80°C over a period of 
10 minutes. 

To confirm the accuracy of genotyping by this 
method, we randomly selected 50 DNA samples and 
subjected them to polymerase chain reaction and re- 
striction fragment length polymorphism analysis or to 
direct DNA sequencing. In each instance, the genotype 
determined by the allele-specific DNA primer-probe 
assay system was identical to that determined by the 
confirmatory methods. 

Statistical Analysis 

Hardy-Wernberg equilibrium was tested separately 
in the patients and controls using ^-tests. To analyze 
interactions between SNPs and between SNPs and 
asthma development risk, the multifactor dimension- 
ality reduction (MDR) approach was used. 11-14 MDR 
was implemented with a cross-validation framework 
(typically, 10-fold, although any number of intervals, 
including 1-fold, can be used). In step 1 of the process, 
the data were divided into a training set (nine-tenths of 
the data) and a testing set (one-tenth of the data). Next, 
a set of "n" genetic and/ or environmental factors was 
selected from the set of all possible factors (i.e., in the 
study n was 1, 2, 3, and 4 for SNPs); each of n individ- 
uals in the training set was grouped according to its 
state at each of the n factors (i.e., for two loci with three 
possible genotypes, there were 9 possible combina- 
tions; for three loci, there were 27, etc.). Then, each 
genotype combination was classified as high or low 
risk depending on the ratio of cases to controls with 
that genotype combination. If the ratio was <1 (i.e., 
there were less cases than controls), the genotype com- 
bination was classified as low risk. If the ratio was >1, 
the genotype combination was classified as high risk. 
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Table 2 SNPs included in the MDR analysis 



Gene 


SNP ID 


Chromosome 


Polymorphism 


Allele 


Minor Allele 

Frequency in Frequency in 
Asthmatic Controls 
Patients 


p Value of HWE 
Test 


CysLTRl 


rs321029 


X 


G>A 


A 


n 4R 


n 46 


<0.0001 


IL-5 receptor 


rs334787 


3 


T>G 


G 


U.IO 


fl ?fl 
U.Z.U 


0.849 


MMP8 


rsl7099451 


11 


OA 


A 


n 1 7 

U. 1/ 


n ?s 


1.000 


B2AR 


rsl042713 


5 


Argl6Gly 


A 


U.JU 


n 4.7 


0.444 


CTLA4 


rs231775 


2 


G>A 


A 


0 38 

\J.OO 


n 4^ 


0.819 


IL-3 


rs40401 


5 


T>C 


C 


U.rtU 


fl 4Q 


0.739 


CRP 


rsl 130864 


1 


C>T 


T 


u.uo 


n ns 

U.UO 


0.629 


CYP2C9*3 




10 


A>C 


C 


0.04 


0.04 


1.000 


CYP3A4*16A 


rsl2721627 


7 


C>G 


G 


u.uz. 


n n? 


1.000 


CYP3A4*1G 




7 


IVS10 + 12G>A 


A 


n 93 


n 99 


0.281 


CYP3A4*11 




7 


T363M 


T 


n nn 


n nn 

u.uu 


1.000 


CYP3A4*18A 


rs28371759 


7 


T>C 


C 


n m 

u.u± 


n m 


1.000 


ADAM33 


rs3918395 (M +1) 


20 


G>T 


T 


U.JLO 


n 14 

U. 14 


0.154 




rs3918396 (SI) 


20 


G>A 


A 


0.00 


0.00 


1.000 




rs2787094 (V4) 


20 


G>C 


C 


0.35 


0.33 


1.000 




rs511898 (F +1) 


20 


C>T 


T 


n 36 


n 49 


0.908 




rs528557 (S2) 


20 


C>G 


G 


n 9R 


n 97 


0.148 




rs44707 (ST +4) 


20 


A>C 


C 




n 46 

U.4U 


1.000 




rs597980 (ST +5) 


20 


C>T 


T 


n 3Q 


n 4n 


0.565 




rs574174 (ST +7) 


20 


C>T 


T 


fl 1 9 

U. -LZ. 


n 1 ^ 


0.800 


CysLTRl 


rs321081 


X 


G>A 


A 


0.44 


0.44 


<0.0001 




rs320995 


X 


T>C 


C 


0.44 


0.44 


<0.0001 




rs321092 


X 


A>G 


G 


0.45 


0.46 


<0.0001 


CysLTR2 


rs41347648 


13 


A>G 


G 


0.00 


0.00 


1.000 


ECP 


rs2073342 


14 


Arg97Thr 


C 


0.23 


0.20 


0.873 


GR 


rs6195 


5 


A>G 


G 


0.00 


0.00 


1.000 




rs6191 


5 


G>T 


T 


0.19 


0.20 


0.292 




rs258750 


5 


T>C 


C 


0.06 


0.09 


0.693 




rs852978 


5 


T>C 


C 


0.04 


0.07 


0.264 


LTA4H 


rs2660845 


12 


G>A 


A 


0.45 


0.41 


0.911 




rs2247570 


12 


A>G 


G 


0.13 


0.10 


0.595 



IL = interleukin; MMP8 = matrix metalloproteinase 8; CTLA4 = cytotoxic 7 '-lymphocyte-associated antigen 4; CRP = 
C-reactive protein; ADAM33 = a disintegrin and metalloproteinase; CysLTRl = cysteinyl leukotriene receptor 1; CysLTRl = 
cysteinyl leukotriene receptor 2; ECP = eosinophil cationic protein; GR = glucocorticoid receptor; LTA4H = leukotriene A4 
hydrolase; HWE = Hardy-Weinberg equilibrium; SNP = single nucleotide polymorphism; MDR = multifactor dimensionality 
reduction; B2AR = j82 adrenergic receptor. 



This collection of high and low risk groups formed the 
MDR model for that combination of factors, and the 
fitness of that model was determined by classification 
error (the number of individuals incorrectly labeled by 
the model). To evaluate the predictive ability of a 
model, its prediction error was calculated based on the 
proportion of individuals in the testing set that were 
labeled incorrectly by that model. Each combination of 
variables specified by the user was evaluated in this 
fashion, and the entire procedure was performed 10 
times, using one-tenth of the data as a testing set and 
nine-tenths as a training set. To evaluate MDR analysis, 
we evaluated balanced accuracy ([sensitivity + speci- 
ficity])/^). 15 The analysis was repeated 10 times after 



reordering individuals, and the averages of cross-vali- 
dation consistency (CVC), training, and test accuracies 
are presented. To decrease random effects due to the 
arranging of data sets, the 10-fold cross-validation was 
repeated 10 times using different random numbers. To 
avoid the effect of linkage disequilibrium (LD) among 
SNPs in the same gene, the SNPs of which the pair-wise 
LD coefficient r 2 exceeded 0.8 were omitted from the 
analysis. For the SNPs that had been selected by MDR 
methods, we compared differences in the allele fre- 
quencies and genotype distribution of each polymor- 
phism between the case and control subjects by using a 
2x2 (allele) or 2 X 3 (encompassing dominant, 
codominant, and recessive genotype models) contin- 
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gency ^-test, with 1 or 2° of freedom. For the afore- 
mentioned analyses, the exact p values were calculated 
using Fisher's exact test. The empiric p values were also 
calculated by permutation tests (10,000 permutations) 
for the allelic analyses. All of the analyses other than 
MDR analyses were performed with PLINK Version 
1.0.6 (MIT Broad Institute, Cambridge, MA). 16 

RESULTS 

All SNPs were in Hardy- Weinberg equilibrium (p > 
0.05) in both the study groups, except for the CysLTRl 
polymorphism, because it was located on the X chro- 
mosome (Table 2). 

MDR Analysis 

Ten repetitions of 10-fold cross-validation resulted in 
a mean CVC of <30% when a model involving a com- 
bination of three or more loci could not detect a signif- 
icant gene model. An SNP, rsl7099451, in MMP8 was 
selected by the single gene locus model; the cross- 
validation mean for this model was 87.0%. The two- 
gene loci model selected two combinations, viz., MMP8 
and ADAM33 SNP rs44707, as well as MMP8 and IL-3 
SNP rs40401, of which the CVC means reached 45.0% 

Using the aforementioned selected SNPs, the single 
gene locus model was found to be superior to the 
two-locus model, based on balanced accuracy and sen- 
sitivity evaluations, when used to predict bronchial 
asthma development in the testing data. The balanced 
accuracy means for these SNPs were 55.7 and 54.5%, 
and the sensitivity means were 68.7 and 56.2%, respec- 
tively. The mean specificity values were 42.7% for the 
single gene locus model and 52.8% for the two-gene 
locus model (Table 3). 

Association Analysis of MMP8 and IL-3 Genotypes 
with Asthma 

For the three SNPs selected by the MDR method, 
association analysis of the MMP8 rsl 7099451 and IL-3 
rs40401 yielded the most significant result using the 
2x2 dominant model (Table 4). The minor alleles 
were A and C in the MMP8 and IL-3 SNPs, respectively 
(Table 1); the risk alleles were the major alleles. 

Association Analysis of MMP8 and IL-3 Alleles 
with Asthma 

MMP8 rsl 7099451 was associated with an increased 
risk of asthma development (odds ratio [OR], 1.64; 
95%CI, 1.07-2.51; p = 0.012; Table 5). The empiric p 
value calculated by permutation testing was 0.014. In 
the same manner, IL-3 rs40401 was associated with an 
increased risk of asthma development (OR, 1.44; 
95%CI, 1.05-1.97; p = 0.025). The empiric p value for 
this association was 0.049. 



Association Analysis of Each MMP8 and IL-3 
Haplotype with Asthma 

According to the gene database, 17 MMP8 rsl7099451 
was the best tagSNP for the LD bin that included 
MMP8, and IL-3 rs40401 was located within an LD bin 
that contained IL-3 (Fig. 1). Each of these SNPs was 
used to construct a haplotype with neighboring mark- 
ers with which it was in LD; in particular, IL-3 rs40401 
was in complete LD with another IL-3 SNP, rs31480. 

Haplotype analysis revealed that distributions of the 
MMP8 rsll225395 and rsl7099451 GC haplotype and 
that of the IL-3 TT haplotype were significantly higher 
in asthma patients in comparison with controls (OR, 
1.52 and 1.43; Table 6). 

DISCUSSION 

In this study, we investigated whether any of 31 
SNPs in genes encoding cytokines or monokines are 
associated with asthma development in a Japanese 
population. In a single locus model, we found that 
MMP8 rsl7099451, and in a two-locus model, the com- 
bination of this rsl7099451 in MMP8 SNP and rs44707 
in ADAM33, as well as this MMP8 SNP and rs40401 
IL-3, were selected as significant asthma development 
risk markers. 

The single locus model of the MMP8 was consis- 
tently selected, and the balanced accuracy and sensi- 
tivity evaluation against the test data were superior 
in the single locus model compared with those in the 
two-locus model. In these analyses, MMP8 rsl7099451 
and IL-3 rs40401 had the most significant effect on 
asthma development according to the 2X2 dominant 
model. It is known that MMP8 is involved in the patho- 
genesis of some diseases, such as cancer, chronic ob- 
structive pulmonary disease, and bronchiectasis. 18-20 
Until recently, MMP8 has been thought to be produced 
only by neutrophils, but it is now known to be ex- 
pressed in some other cell types, including articular 
chondrocytes, fibroblasts, endothelial cells, corneal ep- 
ithelial cells, smooth muscle cells, and cytotropho- 
blasts. 21 " 25 

Using an MMP8 knockout mouse model, Gueders et 
al. revealed an important role for MMP8 in the control 
of neutrophilic and eosinophilic infiltration during al- 
lergen-induced lung inflammation; they also showed 
that the anti-inflammatory effect of MMP8 is partly 
caused by regulation of inflammatory cell apoptosis. 26 
In asthma patients, a significant inverse correlation 
between bronchoalveolar lavage fluid MMP8 levels 
and FEV X and bronchoalveolar lavage fluid-activated 
MMP8 forms and FEVi was detected 27 Inhaled budes- 
onide normalized the MMP8/ tissue inhibitor of 
MMP-1 ratio in asthmatic children via up-regulation of 
tissue inhibitor of MMP-1 production and down-regu- 
lation of MMP8 production by airway macrophages. 28 
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Table 3 Results from MDR analysis 



Locus 
No. 


Predictive Marker 


Average 
Cross-Validation 
Consistency (%) 


Average Balanced 
Testing Accuracy (%) 


Average Testing 
Sensitivity (%) 


Average Testing 
Specificity (%) 


1 


MMP8* 


87 


55.7 


68.7 


42.7 


2 


MMP8* and ADAM33# or 


45 


54.5 


56.2 


52.8 




IL-3§ 










3 


Not determined 










4 


Not determined 











*MMP8 (rsl7099451). 
#ADAM33 (rs44707). 
§UL-3 (rs40401). 

MMP8 = matrix metalloproteinase 8; ADAM33 = a disintegrin and metalloproteinase 33; IL = interleukin. 



Table 4 Association analysis of MMP8 and IL-3 genotypes with asthma 


SNP ID 


Genotype 


Asthmatic 


Controls 




Analyzing Model 








Patients 




Genotypic 


Dominant 


Recessive 










p Value 


p Value 


p Value 


rsl7099451 (MMP8) 


AA 


7 


6 


0.022 


0.007 


0.570 




AC 


55 


51 










CC 


144 


70 








rs40401 (IL-3) 


CC 


35 


30 


0.073 


0.025 


0.156 




CT 


96 


65 










TT 


75 


32 









MMP8 = matrix metalloproteinase 8; IL = interleukin; SNP = single nucleotide polymorphism. 



Table 5 Association analysis of MMP8 and IL-3 alleles with bronchial asthma 



SNP ID 


Allele 


In Asthmatic Patients 


In Controls 


p Value 


OR (95% C) 


rsl7099451 (MMP8) 


A 


69 


63 


0.012 


1.64(1.07,2.51) 




C 


343 


191 






rs40401 (IL-3) 


C 


166 


125 


0.025 


1.44 (1.05, 1.97) 




T 


246 


129 






MMP8 = matrix metalloproteinase 


8; IL = interleukin; SNP = single 


nucleotide polymorphism; OR 


= odds ratio. 



Overall, these data suggest that MMP8 and its activa- 
tion have an important role in the airway destruction, 
healing, remodeling, and treatment response in 
asthma. 

Our results suggested that the minor A allele of 
rs 17099451 in MMP8 may have a suppressing effect on 
asthma development. Conversely, in a person carrying 
the major C allele of this SNP, the promoter activity 
and the ensuing gene expression may be altered, so 
that suppression of inflammation is lost, resulting in an 
increase in airway inflammation and remodeling and 
the development of asthma. 

MMP8 forms part of a cluster of MMP genes, MMP1, 
-3, -7, -10, -12, -13, -20, etc., located on chromosome 
llq22.3. 29 In this study, we found that rsl7099451 



in MMP8 was in LD with another MMP8 SNP, 
rsl 1225395; it is also possible that this SNP is in LD 
with other SNPs that may contribute to asthma devel- 
opment. 

Furthermore, IL-3 is a cytokine produced by active T 
cells, mast cells, and eosinophils, which stimulates the 
proliferation of hematopoietic progenitors. 30 More- 
over, IL-3 mRNA-expressing cells and IL-3 production 
are increased in asthma patients; this increase in IL-3 is 
independent of eosinophil number. 31 Park et al. 32 re- 
ported that one exonic SNP (in exon 1), IL3 + 79T > C 
(Ser27Pro) showed significant association with the risk 
of asthma and atopy. The Pro-allele had dominant, 
protective effects against the development of asthma in 
nonatopic subjects (p = -.002) and also showed signif- 
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Figure 1. According to the gene database, (A) matrix metalloproteinase (MMP) 8 rsl7099451 was the best tag-single nucleotide 
polymorphism (SNP) for the linkage disequilibrium (LD) bin that included MMP8 (red). (B) Interleukin (IL)-3 rs40401 was located within 
an LD bin that contained IL-3. 
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Table 6 Association analysis of MMP8 and IL-3 haplotypes with asthma 



SNP ID Haplotype r 2 Position Reference p Value OR 

Haplotype 



rsll225395 - rsl7099451 G-C 
(MMP8) 

rs31480 - rs40401 (IL-3) T-T 


0.93 Chrll: 102101689-102101709 (A-A, A-C) 
1 Chr5: 131424230-131424376 C-C 


0.020 
0.025 


1.52 
1.43 


MMP8 = matrix metalloproteinase 8; IL 


= interleukin; SNP = single nucleotide polymorphism; OR = 


odds ratio. 





icant association with the risk of atopy in normal con- 
trol subjects (p = 0.007). 32 We also found that rs40401 
in IL-3 had the most significant effect in the 2x2 
dominant model analysis, which suggests that the mi- 
nor C allele may have a suppressing effect on asthma 
development. 

Van Eerdewegh et al. reported that ADAM33 on 
chromosome 20pl3 was linked to asthma and bron- 
chial hyperresponsiveness, using a genomewide scan 
of 460 white families. 9 The ADAM33 gene group ap- 
pears to be related to proliferation, differentiation, or 
migration of fibroblasts, muscle fibroblasts, or smooth 
muscle cells, factors that are related to bronchial hy- 
perresponsiveness or bronchial remodeling. However, 
these proteins do not appear to be involved in the 
respiratory epithelium, which is related to asthmatic im- 
munologic or inflammatory components. 9,33 Although a 
previous genomewide association study in a Japanese 
population did not implicate ADAM33 among the three 
new susceptibility loci for adult asthma, 34 another study 
suggested the involvement of ADAM33 in the develop- 
ment of childhood asthma among the Japanese. 35 

Our finding that the MDR-based two-locus model 
analysis implicated rs44707 in ADAM33 as a signifi- 
cant factor also supports the relation of this gene to 
asthma development. Moreover, our finding that the 
combination of MMP8 and IL-3 polymorphisms and 
MMP8 and ADAM33 polymorphisms were also se- 
lected by two-locus model analysis may suggest that 
several gene polymorphisms may be involved in the 
onset of asthma. 

CONCLUSIONS 

The finding that an MMP8 allele was most strongly 
associated with asthma indicates that metalloprotei- 
nase function is very important for the airflow limita- 
tion process of the disease. Gene variation may be 
important to the differences in protein expression that 
allow the shedding of cytokines or cytokine receptors 
in this disease. 
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